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ABSTRACT 


“= 1. Introduction 


Auditory physiognomies have been used to localise stenosis in numerous vessels. The non-invasive diagnostic 
technique, if achievable, seems a better approach to the detection of cardiovascular anomalies. This was echoed by 
Akay et al. [1], sought by Lees and Kistler [2] and Miller et al. [3] who used the spectral analysis technique for the 
detection of bruits in the carotid artery. Noticeable diastolic murmurs have also been recorded in patients with 
coronary artery stenosis [4, 5, 6]. Interestingly, Wang et al. [4] reported that acoustic quantities acquired from 
sensitive microphones placed on the chest are being used in a procedure for a non-invasive diagnosis of coronary 
artery disease (CAD). Moreover, the spectral content of isolated diastolic heart sounds has been estimated by 
using exclusively developed signal processing techniques. The result indicates a rise in high-frequency 
components in patients with occlusive coronary arteries. Studies showed that turbulent blood flow generates 
sounds in partially occluded arteries in the event of coronary stenosis [7, 8]. In a physiological state, blood flow is 
largely laminar. Typically, the laminar flow may be disturbed and transit to turbulence mediated by high flow 
velocity in which the critical Reynolds number (Re) is exceeded. Thus, the blood flow may well be described as 
being chaotic as it does not flow linearly and smoothly in adjacent layers. It is pertinent to note that transient 


turbulent flow also applies to the physiological state, as may be observed in many postural maneuvres. 


In the cardiovascular system the vibrating structures, distressed regions of turbulent flow, mixing of flows of 
diverse temperatures, and their semblances produce large pressure fluctuations that are transferred through a fluid 
as sound. The turbulent intensities can be relatively large between the stenosis (where prevalent) and the region 
where turbulence has considerably decayed, and the wall of the artery can be subjected to large fluctuating stresses 
imposed by the turbulent flow. Turbulence is marked by the recurrent unexpected changes in the flow pattern 
which transcends the understanding of the instantaneous motions. The instantaneous fluid velocity induces 


convection which is arguably the cardinal process in turbulent flows. Since molecular diffusion does not give an 
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inspiring contribution to spatial transport at a high Reynolds number, convection mediates the transference of 
momentum, chemical species, and enthalpy [9]. Sometimes a statistical description of turbulence which comprises 
a probability distribution for stationary flows to determine the nature of turbulent velocity fluctuations holds well 
for the description of turbulent motion [10]. Monin and Yaglom [11] articulated the velocity structure functions 
(VSF) which are directly related to the probability density function (PDF) of the local dissipation rate. It is 
therefore a pertinent descriptor of the PDF of the dissipation rate which completely describes the statistics of 


turbulent velocity when considerations are on the Eulerian/Lagrangian fluid particle velocity framework. 


The incidence of turbulence is felt at various segments of the cardiovascular system. In the ascending aorta, the 
laminar flow may transit to turbulent in the event of extremely high flow velocity [12] which is marked by 
intra-arterial ejection murmurs [13]. The coronary arterial blood flow is maximum during diastole, and the 
murmurs induced by turbulent blood flow through partially occluded (stenosed) coronary arteries are most audible 
[1]. Fredberg [14] investigated two major descriptors of the fluctuating turbulent pressure at the wall: the 
root-mean-square pressure and the spectral density of pressure. Elsewhere [15] the spectral analysis techniques for 


blood flow velocity were detailed. 


In this work, we sought a turbulent pressure source point; we regard a single point source in the arterial 
compartment as a source whose dimensions are much less than the wavelength of sound generated. It is to be 
accepted, however, that a random point source is an avid approximation since the lengths of vascular vessels may 
be considered ‘short’. Details concerning turbulent fields in short vessels are sparse. The identification of the 
source point would aid in localizing vascular sounds. The hemodynamic pressure fluctuation on the aortic lumen 
boundary was discussed. It is a factor in the turbulent flow that is often implicated in the creation of aortic stenosis 


murmur. 

oe 2. Fully developed pipe flow 

Fully-developed flow describes the case in which the velocity profile is invariant along the direction of flow (in the 
axial direction). Thus, two different axial locations in the pipe manifest the same velocity profile. Laufer [16] 


investigated the structure of turbulence in fully developed pipe flow. Consider a cylindrical pipe analogous to the 


arterial/venous configuration. 

2.1. Perturbation equations in cylindrical coordinates 

Considering a section of the pipe governed by incompressible mean flow, the continuity equation and Reynolds 
equation in cylindrical polar coordinates may be furnished by Reynolds equations of pipe flow [16, 9]. Using the 


basic Navier-Stokes equations the perturbation equation of continuity reads: 


0p _ 0p Ou; 
+u; —+ “=0, 1 
Ot o OX; - Ox, @) 


where p is the fluid density, and u; is the velocity. The momentum equation for the coherent perturbations is of the 


form: 
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where jz is the dynamic viscosity, and jug is the bulk viscosity (see Pierce [17]) which is assumed constant here. The 


Be 2) 
Ox j , 


perturbation Reynolds stress reads: 
iy = (ulu’,) - uu. (3) 


From Carey [18], the quantity ij may be understood as the oscillation of the background Reynolds stresses due to 


the passage of the coherent perturbations. Equations (2, 3) describe the motion of a Newtonian fluid, accompanied 
by the motion of coherent perturbations, for example, sound wave propagation. In the treatment of cardiovascular 
flows, it is frequently desirable to scrutinize the flow and the coherent perturbation quantities separately, basically 


by some splitting technique (see Morfey [19]). The following decomposition ensues: 
W(x,t) =W(x) +W(x,1) +W(x,0), (4) 
where W encodes the mean (time-averaged) contribution, W is the periodic wave (coherent perturbation) driven 


by the pulsatile structure of the flow, and W signifies the turbulent fluctuation. Now, W indicates any one of the 


variables: 2, u;, T, or p, as the case requires. We define W(x,t) as 
a 1 t+é 
W(x,t) = lim — | W(x,t')dt', (5) 
EDO € ot 


and the phase average is 


N 
. 1 B 
(W(x,t)) = Jim dG t +8) (6) 

P &=0 
where f signifies the frequency of the coherent perturbation. 
A fully developed turbulent pipe flow is characterised by a homogeneous stream-wise mean flow with the velocity, 


u =[U,(r), 0, 0], which is locally parallel. So, U = [ u. x (7), 0, 0] is parallel to the x-axis and fluctuates in the radial 


direction r. The perturbation equations in cylindrical coordinates are: 


( Z +5) P + & ging lo =0 (continuity equation) (7) 
Po or r red ox 


Dp 6 2a ~ 2. I~ 
(2+3)i,+ on LB Sh SH) 


ot Po Po Ox arr Orr? OP ax? 
1 wp \( , 10 ity 107%, Ora, (x—momentum equation) (8) 
+v t t t t 
3 w j)\ @xdr réx0O r Ox ax? 


OF ne wi OF. LF )-& 
dr ri a0) ox 
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(r—momentum equation) (9) 


(@-momentum equation) (10) 


Gaia a — ([Z+3)pec9| i{ : ax ), (The energy equation) (11) 
(oa) evn 
where, 
wan, © and F=0,2, and qa ee, gp _ le Oe 
a ic ar ar ar ox 


“© 3. Flow in Stenosed Arteries 


In the course of systole, blood immediately proximal to a stenosis undergoes laminar flow and experiences 
convective acceleration as it passes from the intact portion of the artery through the converging section of the 
stenosis. At a reasonably high Reynolds number, the flow through the diverging section of the stenosis split up 
from the walls owing to its inability to overcome the adversative pressure gradient. A shear layer susceptible to 
fluid-dynamical fluxes is created at the frontier between the separated high-velocity jet and the rather slow-moving 
fluid in the recirculating separation zone. The shear layer created delivers a source from which these instabilities 
derive energy from the mean flow [14, 20, 21]. This energy extraction process ensues at an amply rapid rate that 
the instabilities break down into fully turbulent motion before the end of systole if the jet Reynolds number is 
sufficiently high. As soon as the jet fills the artery, the turbulence begins to evanesce as it could no longer extract 
energy from the mean flow. The turbulent intensities may be reasonably large within the stenosis and the region 
where turbulence has considerably decayed. Consequently, the conducting arterial segment may be exposed to 


enormous fluctuating stresses due to the turbulent flow [14]. 
3.1. Noise (sound) generation 


Noise is a sacrosanct element of turbulence. Acoustic sources which radiate sound develop from velocity and 
pressure fluctuations. Such sound has varying intensities and directionality that are characteristic of the precise 
nature and magnitude of the sources and the acoustic features of the adjacent media. Pressure is a factor of interest 
in the wave equation because ears and the auscultation devices are sensitive to pressure fluctuations. Turbulent 
fluctuating pressures which are related to acoustical pressure are called sound pressure or sound, and intra-arterial 


sounds are fluctuations of pressures that are related to sound. [22]. 
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Acoustic velocity is commonly used as a boundary condition in computations involving solid bodies. It is the usual 
fluid dynamic boundary condition that must be satisfied [17]. Thus, a good relationship is found between pressure 
and velocity. Morfey and Wright [19] formulated bounded domains that enable the application of pressure-based 
options to the fluid density as a wave variable, thereby refining the pioneering work of Lighthill [23] on acoustic 


analogy. 


The expressions for conservation of mass and momentum of the windowed equations of motion for a fluid 


occupying a region, R., are 


é to} to} 
(pu; ) =0, ay (Pui) ay (eit + pi) = 8, (12) 


a) 
(e pote 


0 


In (12) above and afterward, an overbar implies that the quantity or variable is being multiplied by some Heaviside 
functions H and &; it is windowed in space and time; p,; =P; — jd, where P is absolute pressure and P;is the 
compressive stress in the fluid; 6; is the Kronecker delta and B;is an applied body force per unit volume. The 
Heaviside functions H and © are such that: 


Hm in R andon S (13) 


0 in’ (prime denotes complement); 


51) = 1 fort =0 
"0 for <0 > (14) 


In (13), n encodes a local normal coordinate described for points in the neighbourhood of S$ and n=fS is a moving 


closed surface in an appropriate dimension, assumed smooth. \Vf | is single-valued. The following derivatives 


apply from the property of the Heaviside function (see Morfey and Wright [19]): 


oa) 50), 0a 0, OH 5(n), OH sti (15) 
Ot OX, on Ox, 


where 0(.) is the Dirac delta function and n;= On/ Ox, is the unit normal to S. For any field variable € 


0g Oe | eo OE Gs | . (16) 
eee Sa E[ yn 6()E- 7 6) | 


i 


where v; = vn; v is the normal velocity of the surface S directed into R. Apply equation (16) to equation (12) to get 


0 0 (= 

at (p Po) ax, (pu;) = (2-0) 1500) +[ pu, —(p= po) ¥; [nSME (17) 
and 

© (71) + - (pun; + Dj) B; + pu, A 5(0)+| pu,(w; —v;) + py }n SOE. . (18) 


It is plausible to eliminate pu; from (17) and (18) to obtain the second time derivative of the perturbed density 


valid for all (x; £). Accordingly, 


2, 
a c 
_ (p po) - (ep pr) A 501-= [pu 600) +< fT pu, —(p- po); [mde | 
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a _) 6B @ ;——— 19 
Be Lm os») + 4 mon |S ga (emma) ( ) 
The density form of the acoustic analogy reads: 
1 2(———)]_ 2 é a = 
(55-" J (o-»») | 7 gle PH 50-5 Lou AH 60 ]+ [1m 5E] (20) 
é 1B, OT: 
Fy LEvrtsS Dz | - at Maes 


Tj, Ji, and Lon the right of (20) encode: 
Lighthill stress tensor 
T, ; = puu; + Py — C5 (P— Po) 5; » (21) 


Surface mass flux vector 


J = pu; —(P- Po)¥; = (Mi —V;) + Povis (22) 
Surface momentum flux tensor 
Ly = pu; (u;—v;)+ Pij- (23) 


The right-hand side of (20) represents the sources that characterise the physical properties of the acoustic wave. 
The first two members are the combined impulsive mass and momentum that drive flow from its original reference 
state; the second and third members contain Ffowcs Williams & Hawkings (FWH) formulation [22] of surface 
monopoles and dipoles, demonstrated by ©, and subsequently, the volume source terms, with the body force B; and 
the stress tensor 7; represented spatially and temporally by H &. Reethuf [24] advocated that mass-flow 
fluctuations may be associated with monopole sources as compact fluctuating forces, momenta, or pressures in a 
flow field relate with dipole sources; and fluctuating shear stresses relate with quadrupole sources. The FWH 
equation (24) with no initial source terms provides an extrapolative bridge between the simulation domain and the 
acoustic far field. The substitution of the local density p by a new variable p* related to p and making some 


arithmetic operations give the equation for g» ( p'-2 ) /or? in the form (see Reneman et al. [25]): 


O(p--) 2 a a 
oe (0° - 2) 501-5 |p uA ote) |+—[ Fino | 
2 [inane] SS -(p-p") a +a oe (rau, +py) ? (24) 
where ;- and ;. are similar to J; and L; when p’ replaces p and 
ee ee ae (25) 


Dt 


A general acoustic analogy now reads: 


ISSN: 2581-5059 OPEN @ACCESS 159 


& 


PROMOTING RESEARCH & KNOWLEDGE 


S Mediterranean Journal of Basic and Applied Sciences (MJBAS) 
Volume 7, Issue 1, Pages 154-166, January-March 2023 


{ Le va (-" =r) = aC = py) H 610) ]-—[p'u. AH 51) ]}+ [Jn 50] (26) 


Gi rs _ > ap, ) &Tr*. 
-2finowe]+ : 6 wy, +2 B+2 p ey ipa 
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where the penultimate term is got by letting the equation of conservation as 


Be 1 Pi (27) 


and 7° is defined in the same way as T;; when p’ replaces p 


3.2. Vascular sound field from a source point 


For ease of analysis, but with no loss of generality, we conceive of a single point source in the arterial compartment 
as a source whose dimensions are much less than the wavelength of sound generated. The lengths of vascular 
vessels may be considered ‘short’. Since details concerning turbulent fields in short vessels are sparse, the 
approximation of a random point source holds well. Assume xg € (J), /) is a source point in an arterial segment. We 


define the source point as 
Payy (O.%) = Pay (OLO(X=X,) ? (28) 


where Pjur» iS the turbulent pressure, the angle frequency, f the effective length of the source, and x, is the 
location of an equivalent source. The 6 function here has units of L’' (as in Succi [26]). The spatial Fourier 


expansion of the left-hand side of (28) gives 


F ik Xo 


1 ~ 
»(Q,K,) = dx, sur (OD, Xy JE 
PyylOrk.) =~ | A, Pus(@%) eS 


_i 
21 


L ee 
fax... (a, x, Je’ 

4 

where « is the wave number. J;, is measured as the vessel’s origin at the inception of turbulence, which may be 
numerically zero. It is assumed that the turbulent sources are contained completely inside the vessel and are 
sufficiently far from the ends. So, p;,» alongside the convective derivative dp ju;p/dx9 could be assumed zero at the 
inlet and exit. We situate the observation point outside the source region. The acoustic field of the point source 


reads: 


Pay (Oa, JitKL 1 as 1 teietien tee (30) 
vl, = turb 1 = s + 7 R, TR ) bg Ay . 
p(a@,l,) al) aay (i=) 2 e | 


where M is the Mach number, Rp pressure reflection coefficient at pipe exit, /,. The radiated field is 


: ant (31) 
p(@,x) =" P,,,(@,%,)| K-iM — 
; Ox, 


G,(x|x,), 
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where G denotes Green’s function, which furnishes the linear acoustic response of a pipe and G,,(x|xo ) denotes 
Green's function of the pressure field prompted by a harmonic source of unit strength located at x = xo. The 


averaging over all source positions removes all dependence on x,. Thus, from (30) we get 


eye (32) 
(pul = e fi dx |p(@,l)| = [Pua uaa : + : IR [ 
: ie 4 \e(a)| (I-M) (1-my!* 


Likewise, the average acoustic field of a distributed source over N measurements of incident pressure is 


2 12 2 
(p(t) ) = Lp@., (33) 


KR | [Pav @K I Pl OED ps 
eo) | (~My -My 


3.3. Heart murmur 


A heart murmur is propagated as compression and shear (elastic) waves. The heart may be conceived as an elastic 
spring due to its contractibility during systole and restitution during diastole. It is amenable to structural failure 
arising from the turbulent flow. Therefore, flow dampers may be used to mitigate the discomfort, damage, or 
structural failure resulting from the vibrating tissues. Heart motion may be modelled by using the general Hooke’s 
law in conjunction with the Kelvin-Voigt viscoelastic model. Let us discuss, in passing, the Kelvin-Voigt model 
concerning the heart. (Kelvin-)Voigt model is denoted by a purely viscous damper and purely elastic spring 


connected in parallel as shown in Figure 1 below. 


E 


Figure 1. Schematic representation of the Kelvin- Voigt model [27] 


The total stress will be the sum of the stress in viscous damper and purely elastic spring: 
o,, =0,+0,: (34) 


where subscript S indicates the stress-strain in the spring and subscript D indicates the stress-strain in the damper. 


In a Kelvin-Voigt material, stress o, strain ¢ and their time rates of change are governed by an equation of the form: 


o(t)= Bett) 2. (35) 


Where, in (35) above, E is a modulus of elasticity and 7 is the viscosity. 
If some constant stress 6p is suddenly applied to Kelvin-Voigt material, then the deformations would 


approach the deformation for the pure elastic material oo/E with the difference decaying exponentially. We get: 


é(t)= oA ie), (36) 
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where f is time and tr = 7/E encodes the retardation time. The retardation at some time ¢; (when the material is 


freed) obeys the equation: 
E(t >t, =e(t, err. (37) 


3.3.1. Murmur generation in the aorta 


The governing equations of murmur concerning the Kelvin-Voigt model are (see Bailoor et al. [28]): 


~4+A—+*6 + ip J)_Q 
Fle i I mi (38 a, b) 


ou! Op’ 6 du’ Ou’ 
pitt ave (Sit 
‘Ot Ox, Ox, Ox, O%, 


where p’;; and u’; signify local elastic stress and velocity fluctuations, respectively. The adjoining organs, bones, 
blood volume, and other tissue are denoted by a homogenous medium with density p,, Lamé parameters / and x, 
and viscosity v. The propagation and dissipation of compressive (bulk) and shear waves within this homogenous 
viscoelastic medium are governed by equation (38). The tissue-air interface on the epidermis is modelled as a 


traction-free surface and governing boundary condition read: 


pin, =0 


"eu! . 
GEG eg) ee 
Ox,” Ox, Ox, : 


where n; encodes the surface normal. In elucidating the boundary condition on the aorta-surrounding tissue 


(39) 


interface, Bailoor et al. [28] suggested that the viscous shear stress induced by blood flow may be presumed 


negligible relative to normal (pressure) forces. The consequential boundary condition reads: 


(40) 


where P’=P—P_ is the hemodynamic pressure fluctuation on the aorta lumen boundary. For the homogenous 


medium described above, the compression wave speed is given byc =,/K/p , Where K= A+2u/3 is the bulk 
modulus. Similarly, the shear wave speed is given byc =,/ ul p.- 


In the frequency domain, at low frequencies (<1000 Hz), the solution of equation (38) can be approximated by 


using the free space Green’s tensor: 


ik XxX. 
(ir, @) = : 3H (k r)+(6. -3—)0 Gl) 
G,(©,o) 12n(A +2)! tt, (Kr) + (, A ) 
where 
(yy ik qd) xX, (ab) . 
@= H' ON-  (6H (k))+(6, -3-2 9H, (kr) 
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In (41) above: 


k, =o ( A+2n)/ 9. Encodes a longitudinal wavenumber, 


k =olJulp. is the shear wavenumber, 


H,, is the spherical Hankel function of the first kind, 
x;,; 18 the vector from the source to the monitor point, with magnitude r. 


The superposition principle gives the ‘n’ component of surface acceleration at the monitor points on the aortic 


surface as 
a,(@) =2>(-io) G,,, (r,,@)n,,P(@)AA (42) 


where P(o) indicates the hemodynamic pressure on the aorta lumen boundary in the frequency domain from the 
k” element. 
3.3.2. Sound signals and power 


Turbulent sound signals may be estimated in numerous ways (see Akay ef al. [1], Ahlstrém [29]), Qin [30]). 
Among the time-based characteristics are signal intensity and frequency. One may find the extensive work by 
Allen and Mills [31] on them quite invaluable. The continuous analytic signal consists of the original signal 


alongside its Hilbert transform given by [29, 27]: 
Loe 
See) = sti [Par = 50) +i, > (43) 
att 
where H(t) encodes the Hilbert transform. This signal has the following envelope: 


mS mnrnelt) = ¥|8(0| +|8, 0) - (44) 


The time-averaged frequency spectrum can be produced by using the Fourier transform. The windowed short-time 


Fourier transform (STFT) reads: 
FG. K=Zsowe-Se™ + (45) 


where w represents the time window, ¢ is the translation parameter and k the frequency parameter. Since for the 
STFT the frequency resolution decreases as the time resolution increases and vice versa, a better resolution is may 
be achieved by the use of shorter windows for higher frequencies and longer windows for lower frequencies. Thus, 


we employ the discrete wavelet transform (WT) furnished by Allen and Mills [31]: 


WI(C.7) -Soow[ <5) (46) 


where y is a scale parameter. The mother wavelet denoted w(t) may represent any signal that is decomposed into a 
series of dilatations or compressions. With the apparent loss of the link to local frequency by employing the 


wavelets the scaling term supplants the frequency. 
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On a continuous mode in the time domain, the wavelet transform takes the form: 


waar? fs, | lar - (47) 
y 

We have 

Wt) = Jo sue dt = (sO), ¥, .(0))> (48) 


where the kernel function results from the scale parameter y and time translation, C, of the window function. The 
wavelet transform in the frequency domain can as well be shown (see Qin [30]). The wavelet transforms method is 
crucial in converting the infinite trigonometric function basis in Fourier transform to a finite attenuated wavelet 
basis. It serves a dual purpose of furnishing both frequency and time. Wavelet analysis is an evolving trend in 


signal analysis. 


= 4, Conclusions 


Barring physiological circumstances as may be seen in some manoeuvres and other essences, stenosis is a culprit 
in turbulent cardiovascular flows. During systole, the blood immediately proximal to stenosis is known to undergo 
laminar flow. It only undergoes convective acceleration as it passes from the intact portion of the artery through the 
converging section of the stenosis. Turbulence occurs at a reasonably high Reynolds number and the flow through 
the diverging section of the stenosis splits from the walls on account of its inability to overcome the opposing 
pressure gradient. A shear layer predisposed to fluid-dynamical fluxes is created at the frontier between the 
separated high-velocity jet and the relatively slow-moving fluid in the recirculating separation zone. This shear 


layer that is created delivers a source from which these instabilities derive energy from the mean flow. 


An arterial segment that is admitting a turbulent flow may be exposed to enormous fluctuating stresses. The 
vibrating structures, distressed regions of turbulent flow, mixing of flows of diverse temperatures, and their 
semblances induce large pressure fluctuations in the cardio-vascular system that are transferred through a fluid as 
sound. Studies revealed that turbulent blood flow generates sounds in partially occluded arteries in the event of 
coronary stenosis. Such pathological sounds are referred to as vascular bruits/murmurs. On this note, the salient 


details of this work are: 
(a) The expression representing the sources that characterise the physical properties of the acoustic wave. 
(b) The expression for the vascular sound field from a source point. 


(c) The hemodynamic pressure fluctuation on the aorta lumen boundary and the consequential aortic (heart) 


murmur. 
(d) The analysis of acoustic signals by the suitable transform method. 


Interestingly, as may be deciphered from (d) above, acoustic quantities acquired from sensitive microphones 
placed on the chest are being used in a procedure for a non-invasive diagnosis of coronary artery disease. 
Moreover, the spectral content of isolated diastolic heart sounds/murmurs are being assessed by using exclusively 


developed signal processing techniques. The result indicates a rise in high-frequency components in patients with 
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occlusive coronary arteries. We adduce that Fourier and wavelet transform methods hold well for the analysis of 


sound signals emanating from cardiovascular sounds. 
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